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Abstract 

We calculate the Bx parameter within the framework of the 1/Nc ex- 
pansion. We essentially use the technique presented by Bardeen, Buras 
and Gerard but calculate an off-shell Green function in order to disentan- 
gle different contributions. We study this Green function in pure Chiral 
Perturbation Theory (CHPT) first and afterwards in the 1/Nc expansion 
in the presence of an explicit cut-off to determine Bk and the counterterms 
appearing in CHPT. The high energy part is done using the renormaliza- 
tion group. For the low-energy contributions we use both CHPT and an 
Extended Nambu-Jona-Lasinio model. This model has the right prop- 
erties to match with the high energy QCD behaviour. We then study 
explicit chiral symmetry breaking effects by calculating with both mass- 
less and degenerate quarks together with the real case. A detailed analysis 
and comparison with the results found within other approaches is done. 
Consequences for present lattice calculations of this parameter are then 
obtained. As final result we get 0.60 < Bk < 0.80. If = rud = we get 

'X 
K 



0.25 <B^< 0.55. 
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Figure 1: The AS = 2 Box diagram 

1 Introduction 

In the Standard Model (SM), strangeness changing processes in two units (AS* = 
2) can happen via the exchange of two W-bosons as shown in Fig. |I| (the so- 
called box diagram). This type of transitions contributes to the — mixing. 
Experimentally this mixing is determined by the Kl — Ks oscillation length 
proportional to the mass difference. The — K'^ mixing gives rise to the so- 
called "indirect" CP- violation which is usually parametrized by the CP-violating 
parameter e. Then the Ki state consists mostly of a CP-odd state K2 with a 
small mixing of the CP-even state Ki 

Kl-K^ + sK,, (1.1) 

where \e\ ~ 2.3 • 10^^. And the Ks state consists mostly of a CP-even state Ki 
with a small mixing of the CP-odd state K2 

Ks-K^+€K2. (1.2) 

There are also contributions to this mixing that change strangeness in two units 
through two AS* = 1 transitions separated at long distances. They are important 
to determine the mass difference, e is CP-violating and is dominated by box 
diagram contributions. We will concentrate on those. For an excellent recent 
review on kaon CP violation see Ref. [|l|. 

At long distances, once the heaviest particles (top-quark, VT-boson, bottom- 
quark and charm-quark) have been integrated out, the diagram in Fig. |^ is 
described by the effective AS" = 2 Hamiltonian |^ 

'Hes^^x) = ^(m?,m2,M^,VcKM) Gp asifi)"^ Oas=2{x) (1.3) 
where CA5=2(a;) the following AS* = 2 four-quark operator 

0^s=2{x) ^ <(x)L^,(x) (1.4) 
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with 2Llf{x) = s(x)7^ (1 — 75) d{x) and summation over colours is understood. 
The strong couphng constant as{fi) is the one with three active hght-quark 
flavours. The function JF is a known function depending on the Cabibbo - 
Kobayashi - Maskawa (CKM) matrix elements, top- and charm-quark masses, 
W boson mass, and some QCD factors collecting the running of the Wilson co- 
efficients between each threshold appearing in the process of integrating out the 
heaviest particles. For its explicit form see Ref. 0, is the Fermi coupling 

constant. The global Wilson coefficient is dictated by the anomalous dimensions 
of the operator Oas=2{x). This operator gets only multiplicatively renormalized. 

The matrix element of the operator 0/^s=2{x) between two on-shell kaon states 
is usually parametrized in the form of the i?j^-parameter times the vacuum in- 
sertion approximation (VIA) as follows 

{K'\OAs=2ix)\K^) = ^5^(/i)F|m^ (1.5) 

where Fk denotes the — > fi'^u coupling {Fk = 113 MeV in this normalization) 
and rriK is the mass. The /i-scale dependence of Bk reflects the fact that the 
four-quark operator Oas=2 has an anomalous dimension and its matrix element 
depends on the scale where it is deflned. The anomalous dimension is known 
and using the renormalization group leads to the deflnition of the scale invariant 
quantity 

Bk = BK{lx)asifiT+ (1.6) 

with 

at one-loop. Here P^^^ is the flrst coefficient of the QCD beta function. For 
three active light-quark flavours and Nc = 3 we have a+ = —2/9. Of course, 
the physical matrix element {K^\T-C^^'^\K^) is independent of the scale fi. The 



/X dependence of (|1.5| ) is precisely compensated in TC^^^ to produce a scale 
independent result. It is in this sense that B^ can be considered physical. The 
anomalous dimensions and the extensions to the box diagrams needed are known 
to next-to-leading logarithmic order [^]. We will restrict ourselves to leading 
logarithmic order. Only this order makes sense to the next-to-leading order in 
1/Nc (see below) considered here. 

The vacuum insertion approximation was historically the flrst way this partic- 
ular matrix element ( |1.5|) was evaluated 0. Here by deflnition we have Bk{^j) = 1 
at any scale and we can only obtain an order of magnitude estimate. Next this 
matrix element was related to the A/ = 3/2 part of i^' — > vrvr by Donoghue et 
al.0 using SU(3) symmetry and PCAC. This leads to a value 13 k ~ 0.37. It 
was then found that this relation has rather large corrections0 due to SU(3) 
breaking. Then three new analytical approaches appeared, the QCD-Hadronic 
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Duality approach |Q, QCD sum rules using three-point functions0 and the 1/Nc 
{Nc = number of colors) expansion framework in [|l^. Lattice QCD also started 



producing preliminary results around this time. A review of the situation several 
years ago can be found in the proceedings of the Ringberg workshop devoted to 
this subject [p^ . All these approaches have in common that they try to get a 



numerical value for the parameter and study its dependence on the renor- 
malization scale fi. All of these methods have been updated and refined. The 



QCD-Hadronic Duality update can be found in ||T2[, a QCD sum rule calculation 
IS m PI and the l/N^ expansion method has had the vector meson contribu- 



tion calculated in a Vector Meson Dominance (VMD) model[|l^]. A review of 



recent lattice results can be found in HTSf. A full Chiral Perturbation Theory 



(CHPT) approach to the problem is unfortunately not possible. The data on 
kaon non-leptonic decays do not allow to determine all relevant parameters at 
next-to-leading order {0{p'^)) in the non-leptonic chiral Lagrangian|T^, 0. 
A calculation of these parameters within a QCD inspired model can be found in 
M where the determination of the factor is done to 0{p'^). 

The leading order result for Bk in the 1/Nc expansion is well known 

BK{fi)=BK = ^. (1.8) 

This result is model independent. However, to go further in the 1/Nc expansion 
requires some model dependent assumptions. Different low-energy models are 
then used in variants on the 1/Nc method ITU |. An example is the calculation 



done within the QCD-effective action model | 

In this paper we will use a variation on the 1/Nc method. A first simplified 
version of this calculation has appeared in Ref. There we used the Nambu- 
Jona-Lasinio (NJL) model with four-fermion spin-1 couplings set to zero. The 
conclusion, there, was that although good matching between the cut-off scale de- 
pendence from the low-energy contribution with the perturbative QCD scale de- 
pendence was found, that happens in the region where one expects vector mesons 
to be important. We present now a complete version with spin-1 interactions 
and a much more detailed discussion of the procedure. We use a pseudoscalar- 
pseudoscalar AS* = 2 two-point function in the presence of the strong interaction 
and the effective AS* = 2 action from (|1.3| ). The method and the reasons for 
this are explained in Section |^. In Section ^ we calculate this two-point func- 
tion in standard Chiral Perturbation Theory at next-to-leading order in momenta 
{0{p'^)). This we also use to show how the physical Bk factor can be obtained 
from this two-point function, and the additional information we can obtain from 
our method. Here we also point out the effect of including the singlet r/i. In 
the next Section, ^, we do a first calculation of the non-factorizable part using 
CHPT for the couplings. Then we give a short overview of the extended Nambu- 
Jona-Lasinio model and our reasons for using it. The main part of our work, 
the calculation of this two-point function is described in Section ^. The checks 
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we have on the results are discussed next in Section |^ and finally we present our 
numerical results in Section The conclusions from this work are summarized in 
the final Section ^. Some examples of explicit formulas for some of the diagrams 
appearing are given in an appendix. 

2 The Method and Definitions 

We calculate here not directly the i?x-factor but the AS* = 2 two-point function 

= j <fxe"'-''{Q\T (P^'{Q)P'^\x)T^s=2)\Q) (2.1) 

in the presence of strong interactions. We use P'^^{x) = d{x)i'y5s{x), with sum- 
mation over colour understood and 

rA5=2 = -Gf J dSO^s=2{y). (2.2) 

The reason to calculate this two-point function rather than directly the matrix 
element is that we can now perform the calculation fully in the Euclidean region 
so we do not have the problem of imaginary scalar products. This also allows 
us in principle to obtain an estimate of off-shell effects in the matrix elements. 
This will be important in later work to assess the uncertainty when trying to 
extrapolate from K ^ tt decays to K —>■ 27t. This quantity is also very similar 
to what is used in the lattice and QCD sum rule calculations of Bk- 
The AS" = 2 operator in can be rewritten as 

r^s=2 = -Gp J d'yU^\y)LM 

= -G^/^/rf^xi/rf^X2e--(---)L^'^(xi)L^,(x2). (2.3) 

This allows us to consider this operator as being produced at the ly-boson mass 
scale by the exchange of a heavy X AS = 2 boson. So we first replace the effect 
of the box diagram in Fig. |l| by an effective operator of the type ( p.2|) . This 
then, in order to have a physical definition of the cut-off scale, we replace by the 
exchange of the X-boson. This is depicted graphically in Fig. |^. Notice that the 
identification which fermion is which, is unique in the large Nc limit. 

The Feynman diagram at the quark-gluon level at leading order in 1/Nc is 
in Fig. ^a. The dotted regions are a single quark line filled with leading in 
1/Nc gluon exchanges, a planar diagram. At the next-to- leading order in this 
expansion there are two classes of diagrams. One is the same as in Fig. ^ 
but now there is a non-planar contribution or an extra fermion loop inside each 
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(a) 



(b) 



Figure 2: (a) The operator ( |2.2| ). (b) Its realization via exchange of an X-boson. 
The colour indices that are summed over are marked next to the quark lines. 




X 





(a) 



(b) 



Figure 3: The leading, (a), and the non-factorizable next-to-leading in l/Nc, 
(b), corrections to the AS" = 2 two-point function Yi/^s=2{,'f) and thus to Bk- 
The shaded regions are planar QCD diagrams. The crosses are insertions of the 
pseudoscalar current, P'^^{x). 
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shaded region. These we call factorizable l/N^ corrections. The second class is 
the diagram shown in Fig. ^d. Here the shaded region is filled with gluons in 
a planar fashion. It is this nontrivial class that we will compute in this paper. 
The first class can be calculated completely in Chiral Perturbation Theory. They 
provide the corrections needed to obtain the physical values of Fx, ttik and wave 
function renormalization. All off-shell corrections needed here are also purely 
determined by the strong interaction CHPT coefficients. 

Now we would like to give some arguments in favour of the technique we will 
use. The calculation of the hadronic matrix element in Eq. (|1.5| ) involves the 
mastering of strong interactions at all energies between two very different scales, 
namely, the W boson mass and the kaon mass. This is, of course, where the 
complexity of the calculation arises. Both the quark-gluon momenta and the 
X-boson momentum cover this broad range of energies. While we can use the 
asymptotic freedom of QCD to perform a perturbative expansion for energies 
large enough, we do not know yet how to do a QCD calculation below energies 
around a few GeV. The technique we want to use here is essentially the one used 
by Bardeen, Buras and Gerard in Ref . [1^, [l^ in a slightly different notation and 
making emphasis in the low-energy model to describe the strong interactions. 
We want to impose as many as possible QCD relations on this low-energy model 
(Weinberg Sum Rules and similar relations). More comments on which low- 
energy model we will use and why are in Section ^ 

The crucial point in this approach is that in electroweak matrix elements while 
we cannot keep track of the quark-gluon momenta due to confinement, we can 
keep track of the scale of the operator by looking at the X-boson momentum. The 
same role is played by the photon in the vr"*" — 7r° mass difference case 1 3^, ^ 
Looking again at the diagrams in Fig. ^ and equation ( p.3| ), one can convince 
oneself that the scale dependence imposed by the running of the quark-gluon 
momenta at these energies can be identified with the dependence on the QCD 
renormalization scale under some conditions. First, for the diagrams in Fig. ^ 
one can see that by attaching gluon lines to the quark lines that the dependence 
in a X-boson cut-off can be identified up to corrections of order / with the 
dependence on the gluon cut-off for these diagrams. Here denotes a typical 
external momentum coming through the quarks. The requirement of / /x^ to be 
small is also the requirement that the operator product expansion is still valid 
at the scale /x. Otherwise higher dimension operators are needed to be included 
in the QCD running to take care of the effects of external momenta. Of course, 
these would be the only corrections if the gluon propagator had the perturbative 
behaviour at all scales. This we know is not the case and is the reason why we 
have to go to an effective model at scales below the chiral symmetry breaking 
scale A^. The hope is then that this effective model is sufficiently accurate up 
to a scale fi < where both q'^/fi^ is small and the perturbative evolution has 
set in. In that case a matching between the change in the perturbative part and 
the change in the low-energy part with /i should appear for some scale between 
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and A^. 

We will work in the Euclidean domain where all momenta squared are nega- 
tive. Then, the integral in the modulus of the momentum r in ( |2.3| ) is split into 

two parts, 

i-Mw f^J. fMw 

/ d|r| = / d|r| + / d|r| . (2.4) 
JO Jo J^i 

In principle one should then evaluate both parts separately as was done for the 
7r+ — vr'^ mass difference in the above quoted references. Notice that from the 
diagrams for four-quark operators at quark-gluon level with just one-gluon line 
attached (i.e. order ag), can only generate logarithmically divergent terms in a 
cut-off fi of the gluon momentum. One expects that the same behaviour will 
appear from the low-energy part of the integral in Eq. (|2.4|) for some scale 



/i^ between and A^, as discussed before, when the hadronic interactions are 
included to all orders in momenta. Therefore, here, we will do the upper part of 
the integral using the renormalization group (RG) using the identification of the 
scale dependence discussed above. 

An alternative way of looking at this is to assume that at some intermediate 
scale /i we get from the low-energy calculation 



/o (27r 



(2.5) 



In QCD, the operator T /^s=2{^^) obeys the following inhomogeneous renormaliza- 
tion group equation 



drA5=2(/i) 

d/i 



IAS=2 



-7A5=2(/i) rA5=2(/i) 



TT 




r 



A5=2 



1 



1 



(2.6) 



with 7as'=2 the gamma function of 0/\s=2ix) ^ and a^^^ the strong coupling 
constant to one-loop. This change corresponds to doing the integral in ( |2.5| ) from 
fi to fi + d/i. Integrating then equation (|2.6| ) between the scales fi and Mw, one 
gets the full integral of ( ^l3|) to be the same but integrated up to fi and multiplied 
by C{fi), with C{fi) the corresponding Wilson coefficient. Therefore, 



AS=2 



(2vr) 



sdix2) 



(2.7) 



A strict analysis in l/Nc would correspond to set C{fx) = 1 + d{fi)a'^p{fi) and 
evaluate the second term using factorization in leading l/Nc- d{fi) is what the 
one-gluon exchange diagrams would give with a lower cut-off fi. We will, however. 
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use the full one-loop Wilson coefficient C{fi) = (a^^\fi) / a^^\Mw)'^ ^ ■ Equation 
( |2.7| ), is in fact the equivalent of the effective Hamiltonian in Eq. (|1.3|) . The 
definition of the function includes the factor as{Mw)~°'+ from the Wilson 
coefficient C{fi). This permits, then, the identification of B^ifJ')- 

The basic difference here with respect to the — m^o case is two-fold. 
First the matrix element considered here has anomalous dimensions. It makes 
the identification of the scale dependence with the cut-off in the X-boson mo- 
mentum highly nontrivial. We know that there should be an explicit logarithmic 
fi dependence here. In the vr"*" — tt" mass difference case the matching was of 
order which means that the intermediate momentum regime is less impor- 
tant. The second difference is that in the 7r+ — case the mass difference can 
be related to a vacuum matrix element P3[] . Here this is not possible. So while 
in the other case two-point functions were sufficient we now need to calculate 
four-point functions in the strong interactions. The contribution to the — 
matrix element vanishes in the chiral limit. This, together with the fact that the 
typical scale is the kaon mass makes that the effects of a nonzero quark mass are 
essential here 0. Thus one can only expect matching of the dependence in the 
Bk case (if any) for scales where the effects of nonzero quark masses are small, 
i.e. for scales /x larger than . 

Up to now, we did not need to specify the low-energy model for the strong 



interactions. In Ref. |10|, CHPT was used, however the range of applicability of 



CHPT is precisely below where one can expect a reasonable matching with QCD 



(i.e., above the kaon mass). Then, in Ref. vector meson interactions were 
included using a particular VMD model, namely the Hidden Gauge Symmetry 
model. We will calculate the lower part of the integral by using an Extended 
Nambu-Jona-Lasinio (ENJL) cut-off model and also in CHPT as a test of our 
ENJL calculation. Some reasons for why do we believe this model is more suitable 
for this purpose and its advantages versus other choices are in Sect. ^. Here, 
we only want to point out that ENJL permits the control of chiral corrections. 
The chiral limit is not clear in other approaches when implementing VMD, for 
instance. As a matter of fact, a very important point we want to address in this 
analysis, is the effect of the explicit chiral symmetry breaking. In addition, this 
model allows also for an l/N^. expansion and a chiral expansion like the one in 



CHPT (see Ref. [H). 



*In the electromagnetic mass difference for tlie kaons tfie effects are also expected to be large 
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3 Chiral Perturbation Theory Calculation of 



In this Section we study the AS = 2 two-point function IIas=2{(i'^) in Eq- ( |2.1| ) 
in the framework of Chiral Perturbation Theory; i.e we use CHPT to calculate 
the contribution of the strong interactions at low energies. 



3.1 Lowest Order 



At lowest order in the chiral expansion 0{p'^) ||28| the strong interactions between 
the lowest pseudoscalar mesons including external vector, axial- vector, scalar and 
pseudoscalar sources is described by the following effective Lagrangian 



(3.8) 



where denotes the covariant derivative 



i{v^ + a^)U + iU{v^,- a^) , (3.9) 
and U = exp an SU(3) matrix incorporating the octet of pseudoscalar 



mesons 







^(x) 



A 

V=2' 



V 



K- 



K+ 



\ 



TC 



+ 



V2 Vg 



'V6 J 



(3.10) 



In Eq. ( |3.9| ) f^(x) and a^(x) are external 3x3 vector and axial- vector field 
matrices. In Eq. ( p.8|) x = 2-Bo (-^ + s{x) + ip{x)) with s{x) and p{x) external 
scalar and pseudoscalar 3x3 field matrices and A4 the 3x3 flavour matrix 
Ai = diag(m„, md, m^) which collects the light quark masses. The constant Bq 
is related to the vacuum expectation value 



|gg|0)| 



q=u,d,s 



-F^Bo{l + 0{M)) 



(3.11) 



In this normalization, Fq is the chiral limit value corresponding to the pion de- 
cay coupling F-jr ~ 92.5 MeV. In the absence of the U{1)a anomaly (large Nc 
limit) [^, the SU(3) singlet rji field becomes the ninth Goldstone boson which 
is incorporated in the $(a;) fields as 



(3.12) 
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The effective reahzation of the pseudoscalar current P'^^(x) at low-energies 
can be obtained from the divergence of the ds component the axial-vector quark 
current. Then, 

P'^''{x) = V2FqBqK^ + ■■■ (3.13) 

where we have explicitly given the lowest order term. The operator rA5=2 trans- 
forms as a component of a (27l, 1^) tensor under SU(3)l x SU(3)ij chiral rota- 
tions. The realization of this operator in terms of the relevant low-energy degrees 
of freedom is determined uniquely by its symmetry structure. At leading order 
in the 1/Nc expansion, this operator has the well-known factorizable current x 
current structure 

vTsL, = J dS^,{y)C,23{y) (3.14) 

with = C^^^ + C^^^ + ■ ■ ■ in an expansion in external momenta and quark 

masses. The lowest order is = —i^WDf^U. The coupling 6*27 is a scale 
independent quantity modulating the (27/,, Ir) part under the SU(3)i x SU(3)r 
chiral rotations of the Standard Model Lagrangian. In the large Nc limit G27 = 1. 

In fact, at lowest 0{p'^) in CHPT, the effective realization of the rAs=2 oper- 
ator has the factorizable structure in Eq. (|3.14| ) with C^^\ One can now 



perform the calculation of 11^5=2(0'^) at this order. The only diagram contribut- 
ing at this order is (a) in Fig. ^. Then, we obtain 

nA5.2(g^) = - ,f°i,2 g^g27. (3.15) 

Where mo is the chiral limit value corresponding to the kaon mass. This result 
is equivalent to 

{K'\OAs=2ix)\K^) = mlF^g2j{ii) or sg^ = ^6-27- (3.16) 

The coupling constant 5'27(/^) is the overall factor modulating the (27/,, 1r) repre- 
sentation of which the operator Oas'=2(3;) is a component. The relation between 
G27 and 5-27 (/u) is 6*27 = 6-27 (/u)«s (/")"+ • 



3.2 Next-to-Leading Order 

At 0{p^) we have to consider both tree-level 0{p^) counterterms and loops of the 
0{p^) terms. To define these loops one needs to introduce a subtraction point 
z/^. This is the CHPT subtraction scale and is not related to the scale /x used in 
the other sections. See the more extensive discussion in Sect. ^. 

Let us first look at the pseudo-Goldstone boson loops. Here we will present 
both the results in the octet symmetry approximation, i.e. using the U matrix in 
Eq. ( ^.l(j| ) with no 771, and in the nonet symmetry approximation (or strict large 
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Figure 4: Chiral Perturbation Theory contributions to IIas=2{<1^) (I) • (a) Lowest 
order, (b)-(f) Higher order non-factorizable. A circled cross is an insertion of the 
external pseudoscalar current, P''^'^{x), a dot is a strong interaction vertex and 
the square with the dashed lines to the dots represents the AS" = 2 operator, 
rAS=2- The full lines are meson lines. 
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Nc limit); i.e. using the U matrix in Eq. ( p.l2| ). First let us give the result for the 
octet symmetry case. The contributions from pseudo-Goldstone boson loops that 
are factorizable into two diagrams after cutting the propagator of the fictitious 
A^* = 2 X-boson introduced in Sect. H can be reabsorbed in the corresponding 
0{p^) expressions for mo, Bq, and Fq in Eq. ( |3.15| ). They also give wave function 
renormalization. These are in Figs, (a), (b) and (c) in Fig. ^(plus the symmetric 
ones) . 

In addition, there are also non-factorizable contributions from pseudo-Gold/- 
stone boson loops in Figs, (b), (c), (e), and (f) in Fig. where the boxes are 
the rAS'=2 operator at C(p^). The calculation of these loops gives the following 
results for 11^5=2 (f?^)- For the integrals we only give the MS renormalized finite 
part at some scale u. The diagram in Fig. (b) gives 



-^0-^0 



+ 



- ml p2 _ 



d^p 1 

X ^ [(^^ + ml)ml \n{ml/ u^) + 3{q^ + mjjmj^ ln(my z/^) 

For the sum of diagrams (c) and (f) plus the symmetric one, we get 

d^ p2(^p2 _^ 3^2 _ ^2^^ 



,2 BiF^ 



3(g2 



2 ^2 



2 f 
K) 



BlF^ 



(g2 -m\f WFi 3 



-m 



K 



(3.17) 



f J (27r)4 (p2 _ 
(6g2 + m\) In (m^/z/^) + 3q^] ■ (3.18 



Diagram (d) does not occur in pure CHPT. More about this diagram is in the 
Sects. H and 0. Diagram (e) gives 



BlF^ f d^p p^ 



S{q^-mlfJ {2T,Yp^-ml 



BlF^ 



(2 2 N2l fi 2F2 o^K^'^^ / ^'^) (3-19) 

To this result one has to add 0{p'^) counterterms. The counterterms of strong 
origin are factorizable. The corresponding diagrams are in Figs, (a) in Fig. ^ and 
(d) in Fig. ^, plus the symmetric ones. They can be completely calculated from 



the 0{p'^) chiral Lagrangian classified by Gasser and Leutwyler ||30| in terms 
of the Lj, i = 1, ■ ■ ■ 10 and Hj, j = 1,2 couplings. These corrections can be 
reabsorbed in the corresponding 0(p"^) expression for rriK and wave function 
renormalization. In addition there are also O(p^) AS* = 2 counterterms of weak 
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®— • ®—m—® 

(e) (f) 

Figure 5: Chiral Perturbation Theory contributions to nA5=2('?^) (H) . (a)-(c) 
Factorizable loop corrections, (d)-(f) Counterterm contributions. Symbols are as 
in Fig. ^. See text for more explanation. 
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origin. The corresponding diagrams are in Figs, (a) in Fig. ^, (e) and (f) in 
Fig. ^ plus the symmetric ones. The general structure of these operators was 
classified in Refs. [IG, Tf, IS]. Then the 0{p'^) AS = 2 structure of counterterms 



contributing to the AS = 2 two-point function IlAs=2iq ) is 

rS=2 = -G27 Gf j [d, (x^u + f/tx)^3 {y) {x^u + t/tx)^3 (y) 

+ (x^U - [/tx)^3 {y) {x^U - f/tx)^3 {y) 

+ iD^ F^ cS'-iy) {d.x^U - U^D.x)^^ (y) 

+ D,cS''iy){x^u + u^xXi^}jy) 

+ tD,C^:l''iy)[x^U-U^xX^\^iy) 



+ Dr tr {x^U + U^x) {y)cS'{y)C%{y)] . (3.20) 

These O(p^) AS* = 2 counterterms have also a factorizable part which can be 
obtained expanding the large Nc expression in Eq. ( p.l4| ). These can be obtained 
from the corresponding 0{p) and 0{p^) terms of the left quark currents in 
CHPT. Then at large N^, = AL^/F^ and D7 = SL^/F^ and all the other 
counterterms are zero. Here L5 and L4 are two of the O(p^) chiral Lagrangian in 
the strong sector . These factorizable counterterms can once more be absorbed 
in the corresponding expression for Fk- So altogether the factorizable 

counterterms of both strong and weak origin are the needed ones to cancel the 
(CHPT) scale dependence generated by the factorizable pseudo-Goldstone boson 
loops discussed above, so that the factorizable part is completely calculable and 
^ives the result in Eq. ( |3.15|) with mo —>■ rriK, Fq — >• F^ and (m^ + ms)BQ — ^ mi- 
Let us discuss the non-factorizable or next-to- leading 1/Nc order contribu- 
tions. They will give the non-trivial contributions to Bk- Although the calcula- 
tion can be performed in terms of the couplings Di, i = 1, ■ ■ ■ ,7; these couplings 
are unfortunately not known and therefore it is not possible to calculate to 
0{p'^) within CHPT. Meanwhile these couplings are not available experimentally 
one can try to calculate them. To do that one needs dynamical information for 
which some QCD inspired model can be useful. The corresponding calculation 
in an effective action approach can be found in Ref. |^. In the next sections we 
will use the technique and dynamical assumptions explained in Sect, ^to obtain 
in the large limit these Dj-like couplings appearing to all orders in CHPT. 

In the large Nc limit we have that Di, D2, -D3, D4 — AL^/F^, Dq and 
are 0{1/Nc). 



tin fact, this discussion can be extended to all orders in CHPT. 
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Diagram (a) in Fig. || produces two kaon propagators and the result after sub- 
tracting the part that is reabsorbed in the factorizable or large contributions, 
is 



2 BlF^ ^ 



IGDliiy^) (ml -mlf + 16D^(z/2) 



4 



2 f - 8^ j {2ml + -9 + 4 iDliu^) - 4 



(3.21) 

Diagram (e) plus the symmetric one in Fig. ^ produces one kaon propagator 
and the result is 



9 2 ^27 

— 



(3.22) 



32D^(z/2)m^ + 4D5(z/')g2 . 

And diagram (f) in Fig. |^ gives no kaon propagators and the result is 

?,2BlF^G27Dl{u^). (3.23) 



Where we have only written the MS finite part at some scale v in terms of D\{v'^) 
and L[(z/2). 

Summing up all the contributions, namely factorizable and non-factorizable, 
pseudo-Goldstone boson loops and counterterms, we get the following result for 
the full nA5=2(g') to 0{p^) 



nAs=2(g^ 



-2 



{iris + ^dY{<f ~ Tn^--^"^ 



X 



1 + 



K 



a!!^ + b + c 



(3.24) 



With 



A 



-16 [lQTT^F^]Dl(iy'' 



m 



2 \ 2 



m 



+ ln(m^/z/2) 



1 m. 



4 ml- 



3 m: 



4 ml- 



ln(m^ /z/2). 



B = [lQTx'F^]l2\D',{u')-^ 



+ 4(^DK^')+^SK)-4''' 
1 



2 + 



m: 



3 m 



2 \n{ml/u') + ln(m^/z/') + ^ „ 



HmlJu') + 1 



16 [lQ7:'Fi] D;{u') + 



(3.25) 
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Each of these A, B and C terms must be separately scale independent. The 
divergent part of the 0{p^) couplings in Eq. ( p.20|) was also determined in Ref. 
with 



4-4, 



167r2F2 



^ + ^{lE- ln(47r) - 1} 



(3.26) 



and 7i = -1/24, 72 = 0, 73 = 0, 74 = 3, 75 = 0, 75 = -3/2, 77 = 1. The corre- 
sponding divergent part of the 0{p^) strong chiral Lagrangian were determined 
in Bni with 



.d-4, 



167r2 



_l_ + i{,,_,„(4.)-l} 



(3.27) 



The ones we need are r4 = 1/8 and Fs = 3/8. These divergences precisely cancel 
the ones generated by the pseudoscalar meson loops leaving a scale independent 
quantity. In the large counting A, B and C are 0{1). Notice that in the 
limit irts = one gets A = 0. This means that the subtraction constant 
nA5=2(0) vanishes when rris = rud- Here we have used the Gell-Mann-Okubo 
octet symmetry relation 3m^g = Am^j^ — m^. The exact cancellation of the A 
term for equal quark masses (m^ = md) introduces an unknown error in present 
lattice quenched calculations which are done in this approximation. 

From the explicit calculations one can also obtain that the singlet meson, 
rji does not contribute to IlAs=2i(l'^) , and thus to Bk- The relevant flavour de- 



fied 



ss 



which is octet. So only 



gree of freedom in the zero-charge sector is 
SU(3) breaking effects can be important. In order to estimate these we have 
also performed the calculation in the nonet symmetry approximation; i.e. us- 
ing the U matrix in Eq. p.l2| ). This is the strict large limit and is the 
one analogous to the calculation we shall do afterwards using the ENJL as low- 
energy hadronic model. Also, as in the ENJL model calculation, we have taken 
the qq states as basis for the states running in the loops and not the mesonic 
basis of the lowest pseudoscalar mesons, namely vr, K and 77 mesons. This 
will change some of the scale dependence of the logarithms since the dynam- 
ical fields are not the same. In fact, from the calculation using nonet sym- 
metry and qq states we get the same expression as in Eq. ( |3.24| ) but changing 
m^W ln(m2 /z/2) by (2/3) (2 m2^-m2)i(2) \n{{2ml-ml)/iy^) and m^W \n{ml/i^^) 



by 2m^^^^ ln(m^/i^^). The exponent in brackets is for the change in A, the other 
one is for the change in B. We have used here that the dd meson has the same 
mass as the pion while the ss one has a mass 27711- — m^. This produces that 
in this basis of fields the divergences are 71 = —1/8 and the remaining ones the 
same as given above. 

The different chiral logarithms in the octet case versus the nonet case for the 
A and B terms produces a numerical difference, which for a scale u = nip ~ 0.77 
GeV, is 

^ctet = ^nonet + 0.36 (3.28) 
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and 



B. 



octet 



B 



nonet 



+ 0.32 . 



Let us now go to the B^ parameter itself. After reducing the AS* 
point function above we get to 0{p'^) 



B 



(4) 
K 



4 ^ 



1 + 



B + C} 



(3.29) 

2 two- 

(3.30) 



The difference obtained above for the A and B terms from the chiral logarithms 
in the octet case versus the nonet case enlarge b'^^ by 0.09 in the octet case, i.e. 
no Tji. 



3.3 Large Nc expansion 

In this subsection we would like to perform also an l/N^ discussion in the frame- 
work of CHPT without additional dynamical assumptions. As said before in the 
large N^. limit the T^s=2 operator has the factorizable structure in Eq. (|3.14|) . 



We also know that in the large N,, limit G27 = 1. So to all orders in CHPT and 
leading 

Bk = \. (3.31) 

The next-to- leading order 1/Nc corrections to this result are more involved. 
At lowest order CHPT 0{p^) one still has only the factorizable structure at next- 
to-leading 0{1/Nc) when substituting by C^^^ in Eq. ( p.l4| ). This happens to 
all orders in 1/Nc and lowest 0{p'^). Now G27 contains 1/Nc corrections. Then 
to Oip^) and to all orders in the l/iV^ expansion 

= \g,,. (3.32) 

Donoghue et al. |P used SU(3) chiral symmetry to estimate 6*27 from the mea- 
sured —>■ 7r"'"7r° A/ = 3/2 decay rate which is modulated by the same ^27- 
They estimated it to be 

G27 ^ 0.88-|- ~ 0.49 (3.33) 

and 

5^^- 0.37. (3.34) 
These two results above in (|3.31|) and (|3.34|) are quite well established. Now, 



one can go to 0{p'^) and next-to-leading l/N^.. Again one has the factorizable 
contributions that are reabsorbed in the values of F^, mx and in wave function 
renormalization. But at this 0{p'^), there are also non-factorizable l/N^ correc- 
tions. The result at next-to-leading 0{1/Nc) can be obtained from the previous 
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section, 





167r2F2 



{A + B + C} 



(3.35) 



Where now G27 contains 1/Nc corrections and the A, B and C terms have the 
same expression but calculated to 0{1) in the large Nc. 

From the A, B and C terms one can see that there are five structures of coun- 
terterms to be determined. Namely, G27, ^7 - SL^/F^, Dl + Dl- ALU 
and D2 + -D3/4 i To do that we need dynamical information on the strong in- 
teractions. In the next section we will see how they can be determined in the 
approach explained in Sect. |^ using CHPT to next-to- leading order. Then, in 
Sect. H we will use the full ENJL to obtain the low-energy contribution to the 
rA5=2(/^) function. There, the dynamical assumptions are both in the use of the 
ENJL model and in the identification of the cut-off scale yU with the perturbative 
renormalization scale. We are in this case calculating the AS = 2 two-point func- 
tion at next-to-leading order in the 1/Nc expansion and to all orders in CHPT. 
So, we are in fact calculating all the Z)j-like counterterms that appear to all orders 
in CHPT at leading 0{1/N,). 

Let us now sketch how this can be done to O(p^). In the chiral limit there 
are two counterterms 



In the chiral limit and on-shell there is only G27- Then we can obtain 6*27 from 
the chiral limit for = 0. The slope in this case will give the C term. Outside 
the chiral limit there are three more structures to be determined. These can be 
determined due to their different behaviour. The A term produces a pole in 
the BxifJ''^, q^) form factor at = for 7^ m^. From its residue on can solve 
for the Di counterterm. The B term has two possible structures involving D-j and 
respectively, that can be disentangled from the different m^/m]^ dependence 
multiplying them. 

In Sect. § we will see how an analogous discussion can be done to all orders in 
CHPT. There we will first determine the G27 coupling by going to the chiral limit 
and making a numerical fit to to the ratio between the lowest order to the next- 
to-leading 1/Nc order to a polynomial in powers of starting by Then 6*27 
is the term at g^ = 0. Once we have G27, we do both, for the case of degenerate 
quarks = and the real case 7^ m^, another numerical fit to the ratio 
between the lowest order to the next-to-leading l/N^ order to a polynomial in 
powers of g^ starting in g^^. Here, we are determining the A, B, C , . . . terms; 
i.e. the Z)j-like counterterms. 

•'■Notice that the combinations + and _D2 + -D3/4 cannot be disentangled. The operator 
multiphed by can be rewritten in terms of the others using field transformations. 



and 



^ 4 



(3.36) 
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We have one more general result in CHPT. To next-to- leading order in 1/Nc 
but to all orders in CHPT the diagrams that can contribute are still those de- 
picted in Figs. ^ and ^. The vertices now are the ones appearing in the CHPT 
Lagrangian at all orders. There is still no non-analytic dependence on possible 
even with all possible vertices. The result from the previous section that at order 

the dependence on in Ha5=2 is analytic, is thus true to all orders up to 
next-to-leading order in l/N^. 

4 Chiral Perturbation Theory Calculation in 
the 1/Nc Expansion 

In this section we use CHPT to estimate the low-energy contribution to rAs=2(/^) 
in Eq. (|2.5|) . We will only discuss the nonet symmetry case and using the states 
qq as dynamical basis for the states running in the loops as explained in the 
previous section. As said in Section Q at some intermediate energy region we want 
to identify the dependence on a cut-off /i in the X-boson momentum with the 
QCD renormalization scale dependence. Under the conditions explained there, 
one expects it to be plausible for some value below the spontaneous symmetry 
breaking scale. 

This scale /i is thus not related to the scale p of the previous section, even 
though it appears in a similar fashion in the logarithms. If we wanted to extend 
the analysis presented here going beyond the lowest order in CHPT coupling of 
currents to the mesons a similar scale v would appear. This would then cancel 
the dependence on v of the Li in the higher order CHPT Lagrangian. The answer 
then would be v independent. The scale yU is the upper limit of the integral in 
Eq. ( p.5| ) and would still be present. 

We have chosen to route the momentum in the loop integrals r a,s px = + Q 
where px is the X-boson momentum and q is the external momentum. Any other 
routing will induce similar uncertainties of 0{q'^ / ii^). Then, we do these integrals 
in the Euclidean space cutting-off the loop momentum te for |r£;p > /x^ (where 
the subscript E stands for Euclidean). 

We want to emphasize here that this is not a pure CHPT calculation as in 
Section ^. It contains some dynamical assumptions like that we can reproduce 
the QCD renormalization scale dependence in the 1/Nc expansion with a cut-off 
in the X-boson momentum, i.e. that lowest order CHPT is good enough up to a 
scale /i where we can compare with the perturbative part of the calculation. We 
will make then both an expansion in g^/A^, /^^/A^ and quark masses over the 
same scale A^ in the context of an 1/Xc expansion. A,^ ~ 1.2 GeV is the scale 
of spontaneous symmetry breaking. In the previous section the requirement was 
g^/A^ and mj^/A'^ small. 

In addition to being an estimate of the non-leptonic parameters in the sense 
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of Ref. [|T^], CHPT provides a model independent result that will help us to check 
our ENJL calculation for Q'^/A^ and A^^/A^ small enough. 

In this notation of CHPT the QCD quark current P'^'^(x) couples to the ex- 



(^)] 



32- 



ternal source [p(x)]23 and the left current L'lf{x) to - 

At lowest order in the chiral expansion only the factorizable diagram in Fig. 
^ contributes. This contribution is 0{N^) in the 1/Nc expansion. The result is 



■f — mo) 



■9 



(4.1) 



This result when properly reduced is equivalent to 



(W\0^s=2{x)\K'^) 



m, 



o-' • 



(4.2) 



At the next order we have both factorizable 0{N^) and non-factorizable 
0{Nc) contributions. The ones that are factorizable are given in terms of the 
couplings of the 0{p'^) chiral Lagrangian which was classified by Gasser and 
Leutwyler and loops of the 0{p^) chiral Lagrangian. These are completely calcu- 
lable and we do not need any dynamical assumption, they only involve momenta 
of order q"^ and not of order fi"^. These corrections are precisely the ones that 
redefine the BqFq couplings to its O(p^) expressions and are shown in Fig. (a) 
in Fig. ^ and Figs, (a) to (d) plus the symmetric ones in Fig. ^. Those that are 
non-factorizable are shown in Figs, (b) to (f) plus the symmetric ones in Fig. ^ 

Diagram (b) in Fig. ^ gives the following contribution 



1^ d'^TE 

{q^-miyj (2^ 
BIF,' 1 



{rE - Qe) 



+ 



(g2 -m\f WFi 



/i 



Te + rn 
{q^ + m\) fi^ - m^ln 



r| + 2m\ - ml 



' - ml 
ml 



-{q^ + 2m\ - ml) /i^ - [2m 



m„ 



In 



'yu^ + 2ml ~ 



2m^ 



(4.3) 



Keep in mind that the mass of the ss state is 2mj^ 



ml 



The sum of diagrams (c) and (f) plus the symmetric one in Fig. § gives 



X 



m 



K 



BlFi 



3 {q' 



m 



2 ^12 

k) 



^ dVg {rE + qE? [3g| + {te + Qe? + 
{2^Y 



{{te + qEY + rn 



2 f 
K) 



BlF^ 



(g2 - m\Y IQti^F^ 



2 
3 



^ - (3g2 + ml) F.ifx^ q\ ml) + 3g2m^F2(/i^ q\ ml) 



BlFl 



(g2 



m 



2 \2 
K) 



IQtx^F^ 3 



/i 



6g>2 - 3g^ 
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(g2_ 



m 



2 \2 
K) 



IGtt^Fo^ 3 



m 



K 



The functions Fi and F2 are 



A' 1 

ds^A^"' 
2g2 



s + — 
A 



(4.4) 



(4.5) 



using A = y (s + mj^ — g^)^ + 4sg2. These integrals can be performed analytically 
but the result is very cumbersome. The final two lines in Eq. ( |4.4| ) give the 
expressions for = and g^ = 0, respectively. 

Diagram (d) in Fig. ^ gives no contribution at this order. This diagram 
deserves more attention related with the Fierzed terms [0]. It provides in fact 
another non-trivial check. In the soft pion limit, one can relate the VV (AA) part 
of this diagram with a sum rule to some moment of a VV-AA two-point function 
. These issues are discussed in Section ^ 
Diagram (e) in Fig. H gives 
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BlF^ 



(g2 - m\f IQn^F^ 3 



BiFi 



(g2 



IGn^F^ 3 



2 



m 



K 



m|- In 



m 



(4.6) 



K 



The result in the last line above is exact both for g = and for mj^ = 0. 

The quartic dependence in the cut-off cancels between the different diagrams 
as required by chiral symmetry. The 1/Nc expansion and the dynamical assump- 
tions mentioned above have allowed us to make a full calculation of the two-point 
function IIas=2{(i'^) both at next-to- leading 0{Nc) in the l/Nc expansion and 
O(p^) in the chiral expansion. The factorizable loops have a u dependence that 
cancels as explained in the beginning of this section. The integrals in ( |4.3| ), ( ^4.4| ) 
and (|4.6| ) produce both analytical and non-analytical //-scale dependence. The 
cut-off procedure we followed has produced logarithmic dependence in /i for 
proportional to meson masses (see below). This dependence should cancel when 
the chiral expansion for nAs=2(g^) is considered to all orders giving no contribu- 
tion to the running of Bk{^ji) for large The perturbative running is after all 
independent of the quark masses. Lowest order CHPT as used here is expected 
to loose validity before we reach such scales. 

Here, we are actually calculating the combination of counterterms and chiral 
logs in Eqs. ( ^.36| ) and ( ^.25[ ). So in order to extract the values of the coupling 
constants we should compare the full expression for 11^5=2 (g^) in ( p.24| ), with the 
one calculated here. 
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In the chiral case = ma = the calculation of the next-to leading order 
contribution to 11^5=2(0'^) can be written as (assuming jj? > q^) 



(g2 - mlf 



(4.7) 



This result added to the one in (|4.1| ), then properly reduced and compared with 
the Bk definition in Eq. (|1.5| ) gives 



BK{^\q^) = - 



1 



(4. 



As we already knew from the CHPT calculation in Sect. ^ there is no pole for 
this case. Doing now the general 7^ case we can get the 6*27 coupling and 
the A, B and C terms. Then to next-to-leading order in 1/Nc and to order p^, 
i.e. we have neglected the higher order CHPT terms in the last line in Eq. 



G 



27 



A 



+ 



B 



+ 



C 




+ 



(4.9) 



Clearly the matching between the QCD perturbative scale dependence and the 
cut-off /i here in G27, A, B and C is not good enough. Remember they have to be 
separately scale independent. To improve this situation we need to go to higher 
order in CHPT. This we will do in the next sections. Using an ENJL model as a 
good low-energy hadronic model we will calculate to all orders in CHPT in the 
framework of the l/N^ expansion. This O(p^) will serve as a check of our ENJL 
calculation. In fact, the comparison of this result (G27, -B, and C) with the 
ENJL calculation for g^/A^ and /x^/A^ small enough is good. 

Taking the point of view ofRef. m we can compare Eqs. ( [4. 91) and (|3.25| ) 
to obtain 



0, 



22 



;i67r^Fo^] + ^) = 1. (4.10) 

Notice that this assumes that for ji'^ » mj^, this 0{p'^) calculation is enough 
to reproduce exactly the QCD perturbative scale dependence and thus one can 
choose n = p. This assumption, if good anywhere, would be good for scales 
around or below the spontaneous symmetry breaking scale (~ 1.2 GeV). 



Short Description of the ENJL Model and 
its Connection with QCD 



For recent comprehensive reviews on the NJL model and the ENJL model 
[E^, see Refs. 



34, RH]. Here, we only will give a brief introduction and reasons 



why we choose this model. The QCD Lagrangian is given by 



£qcd 

''-QCD 



'^QCD 



Q {h'' {dt, - iVf, - ia^'j^ - iG^) - (M + s - ip'j^)} q ■ (5.1) 



Here summation over colour degrees of freedom is understood and we have used 
the following short-hand notations: q = (u, d, ; is the gluon field in the 
fundamental SU(Nc) (Nc=number of colours) representation; G^j, is the gluon 
field strength tensor in the adjoint SU(Nc) representation; v^, a^, s and p are 
external vector, axial-vector, scalar and pseudoscalar field matrix sources; Ad is 
the quark-mass matrix. 

From lattice QCD numerical simulations, all indications are that in the purely 
gluonic sector there is a mass-gap. Therefore there seems to be a kind of cut-off 
mass in the gluon propagator (see the discussion in Ref. [0). Alternatively one 
can think of integrating out the high-frequency (higher than A-^, a cut-off of the 
order of the spontaneous symmetry breaking scale) gluon and quark degrees of 
freedom and then expand the resulting effective action in terms of local fields. 
We then stop this expansion after the dimension six terms. This leads to the 
following effective action at leading order in the 1/Nc expansion 



-QCD 



with 



and 



S,P 
ENJL 



£enjl 



-QCD 



ENJL 



+ 



i/K 



(A, 



(5.2) 
(5.3) 



V,P 
ENJL 



Sn^Gy (A, 

iV.A! 



(5.4) 
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Here i,j are flavour indices and \l/i?,L = (1/2) (1 ± 75) ^- The couplings Gs 
and Gy are dimensionless and 0{1) in the 1/Nc expansion and summation over 
colours between brackets is understood. The Lagrangian £qcd includes only low- 
frequency modes of quark and gluon fields. The remaining gluon fields can be 
assumed to be fully absorbed in the coefficients of the local quark field operators 
or alternatively also described by vacuum expectation values of gluonic operators. 
So at this level we have two different pictures of this model. One is where we have 
integrated out all the gluonic degrees of freedom and then expanded the resulting 
effective action in a set of local operators keeping only the first nontrivial terms 
in the expansion. In addition to this we can make additional assumptions. If 
we simply assume that these operators are produced by the short-range part of 
the gluon propagator we obtain Gs = 4:Gv = Ncas/ir. The two extra terms 
in ( ^.3|) and ( ^.4|) have however different anomalous dimensions so at the strong 
interaction regime where these should be generated there is no reason to believe 
this relation to be valid. In fact the best fit to the low energy effective chiral 
Lagrangians of 0{p'^) and O(p^) fS^], is for Gs ^ 1.216 ^ Gy ^ 1.263 (Fit 1 



in Ref. fS^I). The third parameter appearing in this picture is also obtained 
from the same fit, with = 1.16 GeV. The light quark masses in Ai are fixed 
then to obtain the physical pion and kaon masses in the poles of the pseudoscalar 
two-point functions [^. Then we get fn = {mu + md)/2 = 3.2 MeV and = 83 
MeV. 

The other picture is the one where we only integrate out the short distance 
part of the gluons and quarks. We then again expand the resulting effective 
action in terms of low-energy gluons and quarks in terms of local operators. Here 
we make the additional assumption that gluons only exists as a perturbation on 
the quarks. The quarks feel only the interaction with background gluons. This 
is worked out by only keeping the vacuum expectation values (VEVs) of gluonic 
operators and not including propagating gluonic interchanges. Most fits are in 
fact better with this gluonic VEV set to zero and when this is not so the results 
are quantitatively very close to the results in that case. Accordingly, we will take 
this gluonic VEV equal to zero in this work. 

This model has the same symmetry structure as the QCD action at leading 
order in 1/Nc ||3^. Notice that the U{1)a problem is absent at this order [^. 



(For explicit symmetry properties under SU(3)l x SU(3)ij of the fields in this 



model see reference ||3^.) The QCD chiral flavour anomaly can also be consis- 
tently reproduced in these kind of models when spin-1 four-fermion interactions 
are included In the chiral limit, this model (for Gs > 1) breaks chiral 



symmetry spontaneously via an expectation value for the scalar quark-antiquark 
one-point function (quark condensate). We use here a cut-off in proper time as 
the regulator. In the mean-field approximation [|l^ we can introduce the vac- 



uum expectation value into the Lagrangian, via an auxiliary field, and then keep 
only the terms quadratic in quark fields. The -C^jl^'^ above are then equiv- 
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(a) (b) 

Figure 6: The graphs contributing to the two point-functions in the large Nc 
hmit. a) The class of all strings of constituent quark loops. The four-fermion 
vertices are either ^Injl -^enjl Eqs. (|5.3|) and ( ^.41) . The crosses at both 
ends are the insertion of the external sources, b) The one-loop case. 



JDI of the form 



alent to a constituent chiral quark-mass term of the type in 

-MQigLU^qn + qBUqi) [0. 

In this model, two-point functions are given by the general graph depicted in 
Fig. ^. How to sum these kind of strings of bubbles of constituent quarks for 
two-point functions, regularization independent relations obtained in this model, 
the extension of the technique from two-point to three-point functions including 
explicit chiral symmetry breaking, discussion of the Weinberg Sum Rules in this 
model, how VMD works in this model and more related phenomenological issues 
are treated in Refs. E^, E^, E3] and reviewed in [E^ . The general conclusion is 



that within its limitations the ENJL-type models do include a reasonable amount 
of the expected physics from QCD, its symmetries, their spontaneous breakdown 
and even some of its short distance information, as for instance the one embodied 
in the Weinberg Sum Rules. This is a very important point and has been one 
more of the reasons why we have chosen this model. Relations like the Weinberg 
Sum Rules are theorems of QCD and should be reproduced by any reasonable 
candidate to describe the low-energy dynamics. In fact, they are essential in the 
convergence of the hadronic contribution to the electromagnetic 7r+ — 7r° mass 
difference [25|. These relations guarantee the good matching between the low- 



energy behaviour and the high-energy one. Models to introduce vector fields like 
the Hidden Gauge Symmetry (HGS) do not always have this good intermediate 
behaviour. This HGS model was used in ref. to estimate vector meson 

contributions to Bk in an 1/Nc expansion. 

The major drawback of the ENJL model is the lack of a confinement mecha- 
nism. Although one can always introduce an ad-hoc confining potential doing the 
job. We will smear the consequences of this drawback by working with internal 

In the 1/Nc expansion. 



\P 



and external momenta always Euclidean with 
we will also only keep singlet color observables. 

We will use the ENJL model as a model to fairly describe in the large 1/Nc 
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Figure 7: The leading 1/Nc contribution to 11^5=2(0'^) in the ENJL model. Sym- 
bols as in Fig. |^ except that a dot is now a ENJL-vertex in Eq. (|5.2|) and the 



full lines are constituent quark-lines. The box is a rAs=2 operator insertion. The 
flavour is mentioned next to the lines. 



limit the strong interactions between the lowest-lying mesons and, if needed, 
external sources. The model we are using is a tree-level loop model with a explicit 
cut-off regularization for one loop parts. This introduces the physical cut-off A^. 
What we mean by a tree level loop model is the following. A general set of external 
fields is connected via a set of one-loop diagrams with three or more legs (vertices) 
and sums over connections of these vertices by a four-fermion interactions or a 
full chain like depicted in Fig. |^. These are also the contributions which are 
leading in l/N^. Going beyond the tree-level approximation is going beyond the 
large A^^c-limit. It is at this level that the hadronic properties of this model have 
been tested. Also to going beyond this one would have to include other operators 
not suppressed at the next-to-leading order in 1/Nc in the ENJL Lagrangian. At 
that level one also encounters the problem of regularizing overlapping divergences 
in the model. For the purpose we want to apply the model here, namely, for 
calculating next-to-leading l/N^. corrections to hadronic matrix elements of four- 
quark operators to consider strings of bubbles is sufficient for the non-factorizable 
part, see Sect. |^. As discussed there the leading non-factorizable part is calculable 
by planar diagrams. This corresponds to the tree-level loop calculation. 



6 The ENJL Calculation 

In this Section we will explain how the calculation was done in the ENJL model. 
In the large Nc limit there is just one kind of diagrams that can contribute. 
These are depicted in Fig. |^. There the fermion lines are constituent quarks 
with mass and Mg We have the product of two two-point functions. 



namely two mixed pseudoscalar-axial- vector Yi^p{q)ad- These two-point functions 
are connected by the exchange of a X-boson between the two left currents. At 
leading order in the 1/Nc expansion (i.e., at 0{Nc) in this case), each one of the 
two two-point functions are sum of all the strings of bubbles or loops, i.e. one, 
two, ■ ■ ■, 00 . This type of diagrams factorizes when the X-boson propagator is 
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Figure 8: The class of three-point diagrams. Symbols as in Fig. |^ the hatched 
areas are a summation over sets of one-loop diagrams as shown in Fig. |^. 
There is also the same type of diagram with the three-point functions tilted and 
the central propagator dd changed to ss. 



cut. In terms of the two-point function [p^, ^ 



(6.1^ 



with A^^{x) = s(x)7^75(i(a;), and Pds{x) = (i(x)z75s(x). Then, at leading order 
in 1/Nc, we have 

1. 



n 



A5=2 



(6.2) 



The two-point function needed here was calculated in the full ENJL model, in 



the chiral limit in Ref. p4| and in Ref. |2^ for non-zero quark masses. In both 
cases to all orders in the external momenta expansion. The result in Eq. (|6.2| ) 
is equivalent to Bk{h) = Bk = 3/4 once ^/is=2{<f') is reduced to the matrix 
element in Eq. ( [1.5|) . 

At next-to- leading order in the 1/Nc expansion, we have two general kind 
of diagrams. The one depicted in Fig. ^ and the one in Fig. ^. These 
are the two possibilities for the tree level loop diagrams and correspond to Fig. 
Pd. In both cases, when we cut the X-boson propagator, we have four external 
legs. Two connected to the left currents and the other two to the pseudoscalar 
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Figure 9: The class of four-point diagrams, symbols as in Fig. ||. 

sources. These are, then, tree-level constituent quark loop four-point functions. 
The diagram in Fig. || can be seen as the product of two one-loop three-point 
functions with a Pds pseudoscalar leg, a left current L^^ leg and the third leg 
amputated. Then the two three-point functions are glued with a propagator, 
this can be ss or dd and with any kind of Dirac structure. The other class 
contains a one-loop four-point function dressed with two legs connecting the 
pseudoscalar sources Pds and another two connecting the left current sources L^^. 
Therefore, we have to calculate in the ENJL model, at leading order in the l/iV^ 
expansion a generalized (P'^'^L^^P'^^L*'^) four-point function. This implies the 
actual calculation of many one-loop four-point functions since the pseudoscalar 
sources and left currents can mix with the other Dirac structures. We have 
also to calculate all possible three- and two-point functions made out of scalar, 
pseudoscalar, vector and axial-vector currents. This, of course, is the major part 
of the work. For a more detailed explanation of the kind of contributions we 
have to consider and some examples of the diagrams in Figs. ^ and ^ see the 
appendix. 

The calculation of these Green functions is done in the ENJL model. The 
following discussion is for non-anomalous Green functions. As emphasized in Sect. 
1^, this is a model of trees of bubbles where one bubble is consistently regularized 
using a cut-off regulator. This cut-off regulator has to preserve the QCD Ward 
identities. We do that by using a proper-time regularization and imposing the 
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QCD Ward identities, i.e. adding the needed counterterms. Although, proper- 
time regularization breaks in principle these Ward identities one can always add 
the counterterms that restore them. Now, we would like to explain, in more detail, 
the way we have done the regularization here. After the Dirac algebra is done, we 
use relations like 2 r • g = ((g + r)^ — M^) — — (r^ — M^) and = (r^ — M^) + 
to reduce the number of propagators. Then, the least divergent integrals 
are just calculated by the standard way of rerouting momenta and reproducing 
the dimensional regularization result when A^ — > cxd. These integrals are the 
ones where all the Lorentz indices are saturated by external momenta indices. 
This only introduces O ^p^/A^j ambiguities but they are of the same order as is 
inherent in the applicability of a non renormalizable model with just these two 
parameters, and A^. Because of the results mentioned in the previous Section 
we do not expect these corrections to be very important. In the rest of the 
integrals, some of the Lorentz indices are carried by g^j^u. These can in principle 
give rise to divergencies that break the Ward identities. In fact the number of 
Ward identities in each case is sufficient to determine these divergent parts. Thus, 
we use all possible Ward identities to determine these kind of integrals. In these 
way we ensure at the same time that our Green functions fulfill the underlying 
QCD Ward identities. In effect only a subset of the Ward identities is needed for 
this. The remainder forms a check on the one-loop calculations. This procedure 
is equivalent to using the heat-kernel expansion. Applying the prescription to 



reproduce the QCD flavour chiral anomaly consistently given in Ref. |23], we do 
not have, here, to add any counterterm to the Feynman diagram calculation. 

Once we have consistently regularized two-, three- and four-point functions 
we close the X-boson propagator by integrating in the loop momentum r up to 
the Euclidean cut-off /i. The routing of the momenta is the one depicted in the 
figures. As explained in Sect. we reroute the external momentum q through 
the X boson momentum px = r + q. The presence of the cut-off in r breaks 
the translational invariance on this momentum and then the two-possible choices 
Px = f + q and Px = f — q give a numerical difference of the order of q^ j \J? ■ 



7 The Vector- Vector Toy Effective Lagrangian 



Donoghue and Golowich |^ suggested to look at a toy effective Lagrangian that 
is a four-quark operator of the type vector current times vector current. The point 
is that the leading contribution to this type of effective Lagrangian is calculable 
in terms of measurable quantities and it may give some feeling on the low-energy 
behaviour of current times current effective Lagrangians like the Standard Model 
one. However, the fact that the low-energy behaviour of this vector current x 
vector current is controlled by leading order in the 'VjNc expansion measured 
spectral functions and thus reliably calculable makes, at the same time, its low- 
energy behaviour quite different from the left current x left current effective 
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Lagrangian one. In this case the low-energy behaviour is given by next-to-leading 
order at large N^. 

The toy AS* = 1 hamiltonian proposed by [RTI] is 



n 



AS=l 
VV 



x,Mw)TlVr^,,iO)V:, 



i5 



X 



(7.1) 



Where g2 is the 811(2)^ coupling and Gp = {92"^^) / i^^w) ■ "^^^ currents are 
V^^^jj{x) = q{x) [(Aq + iXb)/2] 7^g(a;) with the following normalization tr(AaAfc) = 
26ab, and Aa are the Gell-Mann SU(3) flavour matrices. Then, they wrote down 
a sum rule relating, in the chiral limit, the amplitudes of the K'^ tt~ and 
—>■ 7r° transitions to measured vector pv{s) and axial- vector Pa{s) spectral 
functions. The reduced amplitudes for these transitions (see Ref. ^l|] for details) 
are proportional to 



A 



ds s In 



Pv{s)pa{s) 







IE 



(7.2) 



This sum rule is equivalent to the one obtained by Dass et al. p5| for the 
7r+ — 7r° electromagnetic mass difference. 



327r2F2 



3e2 



^ [ml^-ml,^ = dssln(^— j pv{s)pa{s) 



(7.3) 



Notice, that in the amplitude A the scale Mw cannot be sent to infinity. In 
the electromagnetic mass difference, the scale A disappears due to the first and 
second Weinberg Sum Rules. 

This vector x vector effective toy model in Eq. ( |7. 1| ) is related, in the soft 
pion limit, to the vector (VV) part of diagram (d) in Fig. ^ discussed in Sect. H 
through a SU(3) chiral rotation. Diagram (d) is for the — > transition. It 
is then one of the subdiagrams we have calculated for the '^As=2{<f) in the ENJL 
model. The amplitude in Eq. ( |7.2| ) can also be written in the Landau gauge for 
the W boson as 



(7.4) 



Where nl/'*(s) is defined as 



with V;^(x) 



I d'xe^'^-{^\T{V!^max))\^) 



= 6, 



ab 



{q'q' - q'g'") Il\^\q') + q'q^'IlPiq' 



q{x) 



(7.5) 



7'^g(x). Il^^\s) is defined analogously, 
le axial-vector (AA) part of diagram (d) can be anal- 



In the soft pion limit, t 
ogously related via a sum rule to the same amplitude in Eq. (|7.4|) with the 
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opposite sign. Then at leading order (soft pion limit) the left current x left cur- 
rent effective Lagrangian does not contribute to this diagram and the result is 



Following Ref. |43| we can split the integral above into two regions separated 



by a cut-off /x (0 < /x < Mw)- For very small one can use lowest order CHPT 

ny)(s) - (s) ^ 2^ - 8L10 (7.6) 

where Lio is one of the 0{p^) couplings of the strong chiral Lagrangian |^. One 
can already see that the lowest part of the integral diverges quartically. For scales 
H very large one can use QCD ||2^ 



Then the upper of the integral diverges logarithmically. One can use some kind 
of effective low-energy model to improve the lowest order CHPT behaviour in Eq. 
(|7.6|) like it was done in Refs. |2^, ^ for the electromagnetic pion mass difference. 



There the ENJL model and the constituent quark model, respectively, were used 
obtaining a better matching with QCD. However, the divergences here are quartic 
instead of quadratic as there making it more difficult. In fact, as noticed in Ref. 



31[] the largest contribution to the sum rule above is for the range of energies 
between a few GeVs and 10 GeV. Unfortunately, in this region the constituent 
quark model or the ENJL model cannot help since they are intended for energies 
below or around the symmetry scale breaking ~ 1.2 GeV. The large mismatch 
between the scale dependence of the VV (AA) part of diagram (d) at low and high 
energies reflects the large anomalous dimensions of the vector- vector (axial- vector 
- axial-vector) four-quark effective operators. This points also in the direction of 
getting large effective couplings for vector - vector (axial-vector - axial-vector) 
four-quark operators at low energy. 

Nevertheless, we can still use the fact that the cancellation between the VV 
and AA part of diagram (d) in Fig. ^ is exact in the soft pion limit (i.e., = 

= 0) due to chiral symmetry alone. This cancellation will then be a check of 
chiral symmetry for our ENJL calculation!. The cancellation indeed happens to 
all the values of fi calculated in this work with an accuracy better than 1 %. This 
is one more non-trivial check to add to the ones discussed in previous Sections. 

In our numerical results we can also separate out the result for the VV (AA) 
part and relate the result obtained from our full calculation for = mj^ = to 
the result obtained from 



dss' 



n?^(.)-H«( 



(7.^ 

ENJL ^ 



§This cancellation between VV and AA parts is the same that occurs in lattice numerical 
simulations. The cancellation seen there is thus a consequence of chiral symmetry. 
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in the chiral limit (m^ = 0). This provides an exphcit check of the PCAC relation 
in the ENJL model and a very non-trivial check on our full calculation. 



8 Results 

In this Section we are going to discuss the results we get. As in Ref. 



we 



have studied three cases, namely, the chiral case where we set the quark masses 
to zero; the case with SU(3) explicit symmetry breaking 7^ rrid and the case 
with non-zero quark masses but with rris = irtd. This is the first time that the 
chiral case is discussed separately and to all orders in momenta. This permits us 
to obtain directly the coupling G27 introduced in Sect. ^. The remaining cases 
are important to assess the size of the corrections due to the presence of non-zero 
kaon and pion masses. The case = 7^ is done because present lattice 



calculations are done in that limit. In fact, it has been noticed recently that 
quenched lattice data, which are used at present to predict this type of matrix 
elements, cannot be used to extrapolate to physical light quark masses. This 
leaves the error due to the use of degenerate quarks and quenched QCD difficult 
to estimate when extrapolating to the physical Bk- 

All the values of input parameters used in the ENJL model are given in Sect. 
1^. The physical parameters, which are relevant for this calculation, obtained with 
these values are Fq = 89 MeV, m^^ = 495 MeV and = 135 MeV. 

The procedure we followed to analyze the numerical results is also the one in 
Ref. We fit the ratio between the l/Nc corrections and the leading result 



for a fixed scale fi to a/ q"^ + b + cq^ + dq'^ which always give very good fits. The 
a term is the lowest one allowed by chiral symmetry. We also allowed for one 
more term here {d term) to increase the accuracy of the extrapolation. Here a, 
b, c and d are /i dependent for the BK{fi^,q^) form factor. Had we calculated 
just to these a, b and c terms would correspond to the A, B and C terms 

introduced in Sect. |^. Notice that at the order in 1/N(. we are working CHPT 
predicts that the q^ dependence of Has'=2(q'^) is analytic to all orders in CHPT, 
see Subsect. ^3. So we expect to fit our numbers with this form. 

The values of the q^ external momenta used to make the fit are in the Eu- 
clidean region and below the constituent quark production threshold. As said pre- 
viously, this is done to smear the possible consequences of using a non-confining 
model. Once we have these coefficients we extrapolate the Bk form factor to 
q^ = m\ to obtain the physical Bk- For calculating 13 k we use 0+ = —2/9 and 
a^sKMr) = 0.33 ± 0.03 pl which corresponds to A^ = (215 ± 40) MeV to one 



MS 

loop. When the two-loop MS running for as{n) is done, this value is in agreement 
with the LEP value for ^^(M^)^. Since the main source of error in our calcula- 



^If the lower value ai^^(lGeV) = 0.336 ± 0.011, recently obtained in Ref. El from the T 
system, is used then Ai^L. = (125 ± 15) MeV the values for Bk are slightly larger and less 
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(GeV) 






5-(/i) 




5™ 






0.3 


0.60 


0.51(6) 


0.75 


-0.11 


0.64(8) 


0.74 


0.63(8) 


0.4 


0.49 


0.48(2) 


0.72 


-1.6 


0.70(3) 


0.72 


0.70(3) 


0.5 


0.35 


0.36(2) 


0.67 


-4.1 


0.70(3) 


0.66 


0.69(3) 


0.6 


0.18 


0.20(1) 


0.57 


-7.5 


0.62(2) 


0.56 


0.61(2) 


0.7 


-0.05 


-0.06 


0.46 


-12. 


0.52(1) 


0.42 


0.47(1) 


0.8 


<0 




0.30 


-16. 


0.35(1) 


0.23 


0.27(1) 


0.9 


<0 




0.11 


-21. 


0.13(1) 


0.02 


0.02 


1.1 


<0 




<0 






<0 





Table 1: Results for Bk and Bk- The uncertainty due to as{Mr) = 0.33 ± 0.03 
is indicated in brackets. 

tion is of hadronic origin we prefer to give also the running B^i^Jt) parameter. In 
addition, this information can be used with any more accurate vale of a'^p one 
can get in the future. The errors in Table |^ are only from the uncertainty in a^. 

Let us start discussing the chiral case. The corrections we get for the chiral 
case are very large and negative. Unfortunately not very much else can be said 
due to the lack of stability. We can only give a range of values for the value of Bk- 
For the lower bound, it is clear that below 0.3 GeV the QCD scale dependence 
is not trustable. For the upper bound one should remain roughly below the two 
constituent quark threshold to be safe of possible effects due to the constituent 
quark production. Therefore we propose for the chiral limit the following range 

0.30 <Bji< 0.50. (8.1) 

This corresponds to a very narrow window in energy [(0.3 ~ 0.5) GeV]. As said 
in Sect. |]for the chiral case (m^ = = 0) we are actually determining the G27 
coupling. The values of G27 corresponding to the range above are 

0.40 < G27 < 0.67. (8.2) 

Notice that this value is compatible with the one obtained in Eq. ( ^.33] ) from the 
TT^TT^ AI = 3/2 decay rate. Moving away from = we obtain a value 
for the slope C or for D2 + -D3/4. 

Let us go to the equal mass case {rris = rrid). The results for this case are in 
the 7th and 8th columns. Here, we also get that the corrections to the leading 
1/Nc behaviour are negative. The fits done are good and are compatible with 
the absence of the a term. This was also the case in the chiral calculation of 
Sect. and provides another numerical check on our results. We observe also 

stability is obtained. 
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that corrections due to quark masses are positive and tend to bring the value of 
Bx to its large result. Again, we have not a very good matching with the 
leading QCD logarithmic corrections. Although, is somewhat better than the 
one we got for the chiral case. Though there is a narrow maximum around fi ~ 
(0.4 ~ 0.5) GeV, this is just an artifact produced by the perturbative running of 
as{fi) for values as low as yU ~ 0.3 GeV. Thus, we prefer also not to give a central 
value for which will be misleading. Then, for the lower bound we have the 
same argument as for the chiral case. For the upper bound we can enlarge a bit 
the window with respect to the chiral limit due to the presence of nonzero quark 
masses that makes the stability broader. Then we propose for this equal mass 
case the following range 

0.55 < 5^? < 0.70. (8.3) 

Notice that in the equal mass case there is no shift needed due to the difference in 
octet or nonet chiral logarithms. From the fit we can determine one combination 
of Dj and + D^, namely, the b term. From the same fit we also get a value for 
C from the c term. The difference between this value for C and the one obtained 
in the chiral limit indicates the size of the corrections. 

Let us now comment on the more realistic case, (m^ ^ rrid). This is the one in 
columns 4th, 5th and 6th. In the 5th column we present the -Bx(/^, 9^) form factor 
for = —0.001 GeV^ showing the presence of the pole term discussed in Sects. 
1^ and §;i.e. a 7^ 0. For small values of this pole is the overriding numerical 
contribution. Its contribution to B^ is of the right size as expected for a higher 
order CHPT contribution. The value of this pole term allows us to determine the 
value of Di. Again the matching with the leading QCD logarithmic corrections 
is not very good although better than for the other two cases. Clearly, in this 
approach, the presence of masses stabilizes the Bk parameter. As happened for 
the Gv = case in Ref. [^], but now more pronounced, the a term dominates the 
Bx{fJ',q^) form factor for small values of even though is just a small fraction 
when extrapolated to the physical Bk- From columns 6th and 8th one can 
conclude that having degenerate quark masses does not have much influence on 
the final values obtained for B^- This is something that can be used for present 
lattice data, remember that they are still done in the degenerate quark limit. So, 
from this more realistic case and with the complete ENJL, improving previous 
similar determinations [p!0| , |14[, we will give what we consider the result for B^ 
in the large expansion approach. We also prefer not to give a central value 
for the same reasons as given in the previous case. In fact, we believe that also 
from the results in iTDI, \IM one cannot infer a central value and it would be more 



realistic there to give a range of values. Therefore, we propose from the results 
obtained in this work the following range of values 

0.55 < < 0.70 . (8.4) 

These numbers can also be used to determine a second combination of and 
-D4 + D3. So we have a handle on all separate coefficients that contribute to Bk- 
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Let us give now the explicit values for the DiS counterterms obtained from 
this calculation. As said above they can be obtained from the a, b and c terms 
of the fits we have done for each value of the scale fi. They should be /i-scale 
independent if the matching with the perturbative QCD running was perfect. As 
explained before, we beheve that a safe range of scales to make predictions from 
this ENJL calculation is for in between 0.3 GeV and 0.6 GeV. Then for this 
range of and for the chiral logs scale i/ = mp ~ 0.77 GeV we get 



-5.1 ■ IQ-'^ < 



-2.4 • 10"^ < 
0.3 < 

-0.27 < 



Dl < -3.4- 10"^ 
4 



di + d;- 4-^j < 0.6, 

^7 - 8§1 < -0.24. 



(8.5) 



Here we have used the ^ data. To disentangle the and D4 + 
terms, we have also used the = data and the value of G27 obtained from 
the chiral case data. Then these terms will have higher order {0{p^)) chiral 
symmetry errors proportional to quark masses. In addition, we know that the 
D2 + D3/A counterterm can also be obtained from the chiral limit (m^ = = 0) 
data, then we get 



10-2 < 



^ 4 



< 7 • 10- 



(8.6) 



The difference from the previous determination is huge (one order of magnitude). 
This again indicates that actually explicit chiral symmetry breaking corrections 
are very large. All these determinations have, of course, order chiral corrections 
since we are just using CHPT to order to get the A's. 

The chiral symmetry breaking due to the presence of quark masses can then 
be quantified to be 

1.60 ±0.35 (8.7) 



B 



K 



in nonet symmetry and 



1.8 ±0.4 



in octet symmetry. This chiral symmetry breaking is actually very large. However 
it is of the same order as others found in the same system, for instance we have 



— ~ 1 8 

Tn2 ~ ' 



(8.9) 
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which also appears in the matrix element of the Oas=2{x) operator in Eq. (|Ll 
In fact our results point towards an understanding of the discrepancies between 
previous calculations of the Bk parameter and, in particular, the small values 
obtained in Ref. 0, the QCD-Hadronic Duality estimate H, |l^ and the QCD- 
effective action approach [^, although the error bars in this last one are quite 
large. The lattice results [1^ and the next-to- leading l/N^ result of [1^ give 



larger values and are contained within the results of the present work. The QCD 



Sum Rules estimations |]T3| have much larger error bars due to the uncontrolled 
scale dependence in their calculation of Bk{ij)- Within errors this QCD Sum 
Rules estimation is compatible with both groups of results above. The effects 
due to explicit quark masses are different in all these different calculations and 
the uncertainty due to these effects could not be estimated in any of them. The 
difference between the results for of these calculations is of the size of the 
explicit chiral symmetry breaking effects we have obtained here. 

In general, our results were obtained with a ninth pseudo-Goldstone Boson. 
This is the correct thing to do in a strict l/N^. calculation. One needs to extend 
the ENJL model to lift the r]' to a higher mass. In the present work, if we assume 
that, as in the strong-semileptonic sector, the coefficients are well determined 
by a leading 1/Nc calculation and that the main effect of octet-nonet symmetry 
breaking is in the loop calculation, we can then take the estimate given in Sect. 
1^ for this effect, 

Bj^ = (Bk) + 0.09 . (8.10) 

V / nonet 

To obtain our final range for the Bk parameter we will also add an educated 
guess of the corrections of this next-to-leading in l/N^ calculation. Then 

the range we propose for the chiral case is 

0.25 < 5^ < 0.55 (8.11) 

and for the real uis ^ rrid case 

0.60 < 5k < 0.80. (8.12) 



9 Conclusions 

In the present work we have calculated the Bk parameter defined in Eq. (|1.5| ) 
in the 1/Nc expansion. We have essentially used the technique in Ref. ||10| with 
more emphasis in the low-energy contributions and systematic uncertainties. We 
therefore extended their method to Green functions rather than on-shell ampli- 
tudes. The use of this Green function, nAs=2 allowed us to study several issues 
involved in the calculation of non-leptonic matrix elements. 

We have performed a complete CHPT calculation of the AS* = 2 11^,5=2 (?^) 
two-point function in Eq. ( p.l| ), and therefore of Bk to 0{p^). From this calcula- 
tion we conclude that there is an unknown counterterm whose contribution to the 
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Bk parameter vanishes in the hmit = m^. This is a source of uncontrolled 
error for present lattice simulations which are done in this limit and therefore 
making the estimation of the error for the extrapolation from this case to the 
real case unreliable unless the exact rus 7^ ma is done. In addition the present 
quenched data has another source of uncontrolled error when extrapolating to 
real quark masses |^ . 



Then we, using an explicit cut-off fi for the fictitious AS = 2 X-boson in- 
troduced in Sect. ^, have used lowest order CHPT first and the ENJL model 
afterwards, to compute the low-energy hadronic contributions to this two-point 
function at next-to-leading 0{1/Nc). We have then studied the type of informa- 
tion one can get from this kind of calculations for the counterterms appearing in a 
pure CHPT like the one in Sect. |^. We studied three cases, namely the chiral case 
'^s = "^d = 0, degenerate strange and down quarks and the real rUs 7^ md case. 
The chiral limit is not easy to obtain in other popular non-perturbative methods 
like lattice simulations, QCD sum rules or QCD-Hadronic Duality. However, it 
provides very interesting information since, in this limit, the 1/Nc correction to 
the Bk parameter is correlated with the 1/Nc correction to the octet AI = 1/2 



dominating gg coupling |^ . We have shown that this two-point function in fact 
allows to determine all free parameters in CHPT to order needed for a study 
of Bk. 

In general we obtain somewhat less stability, here, in the complete ENJL 
model than in Gy = case studied in Ref. |2l[]. This is just telling us that, 
as said there, the model with Gy = had stability in a region where the spin- 
1 four-fermion interactions are important and cannot be neglected. The study 
done for Gy = was to all orders in CHPT and now we have also the complete 
model also to all orders in CHPT. This situation is not quite the same as the one 



in Ref. |T0|. There, first lowest order CHPT was used and then vector mesons 



added in VMD model to enlarge the range of energies where the calculation was 
reliable. There, as expected then, more stability (enlarged region) was obtained 
when vector mesons were added. Although the results of our calculation have 
slightly less stability than those in Ref. [0, we have, as said in Sect. |], a model 
for including the vector and axial-vector mesons that matches with the QCD 
high energy behaviour. Namely, our ENJL model possesses the first and second 
Weinberg Sum Rules. Then we believe that the results we get are more realistic 
than the ones in Ref. [1^. This lack of good stability (matching) also tells us 
that, unfortunately, this ENJL model fails to reproduce the QCD perturbative 
scaling. We have forced the model to higher scales to see its behaviour compared 
with perturbative QCD. The cause of the failure can be traced in the treatment 
of the gluon propagator. We have taken it to be local below some cut-off scale 
~ 1.2 GeV. This seems a too drastic reduction for energies above ~ 0.6 GeV. A 
more sophisticated model allowing for a momentum dependent gluon propagator 
might help. Nevertheless, we believe that the qualitative conclusions regarding 
the Bk parameter and the very conservative numerical ranges we propose, and 
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in particular, for IBk and explicit cliiral symmetry breaking obtained here, are 
correct. 

The chiral corrections due to quark masses are large and positive. As remarked 
in |]12[, in the QCD-Hadronic Duality approach there are operators of higher 
dimensions that contribute to Bk- These were partially taken into account by 
using intermediate hadronic states with the physical masses and widths. To do 
better than that would require to derive the effective QCD Lagrangian to all 
orders in the expansion in quark masses which is clearly beyond reach now. We 
have instead used an ENJL cut-off to model the low-energy hadronic interactions 
with the advantage that the effective Lagrangian can be derived to all orders in 
CHPT and therefore give an estimate of those higher dimensional operators. Our 
result in Eq. (|8.8|) gives some hint of why the QCD-Hadronic Duality approach 



in Refs. P, 1^] gives lower values than the lattice results or the l/N^ calculations 



here and in Refs. |T0, 113. In fact the QCD-Hadronic Duality approach gives 
roughly the ones we get for the chiral case, these also coincided with the PCAC 
in the chiral limit determination of ||^. 

With the dynamical assumptions presented in Sect. 0, we have also seen that 
the results for the equal mass case and the real case are not very different. This 
conclusion is very important for the present lattice data since they cannot be 
extrapolated to the physical masses for the reasons explained above. They still 
contain the error due to the use of quenched QCD in the simulations. 

We have also included in our final result the shift due to the nonet approx- 
imation inherent in the 1/Nf. limit as well as an educated guess for the l/N^ 
corrections. Our final result for the parameter is: 

0.60 <]3k< 0.80 (9.1) 

Implications of this for CP-violation phenomenology can be found in several 
references, e.g. @]. 



Acknowledgements 

We thank Eduardo de Rafael for discussions. This work was partially sup- 
ported by NorFA grant 93.15.078/00. JP thanks the Leon Rosenfeld foundation 
(K0benhavns Universitet) for support and CICYT(Spain) for partial support un- 
der Grant Nr. AEN93-0234. 



Explicit expressions 

Here we give some explicit expressions and sketch how the four- and three-point- 
like diagrams are calculated. We will then give an example for one three-point-like 
diagram and another for four-point-like diagrams corresponding to the Figures ^ 
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and ^. The notation for the four-fermion couphngs used here with respect to the 
one in Section || is 



9s = 



9v = 



(A.l) 



The AS = 2 two-point function IIas=2{(i'^) we have to calculate (see Sect. 
1^) involves as intermediate stage the calculation of the generalized four-point 
functions {Pds{Q)L^'^^{x)Pds{y)Lf{z)), with 2L^'{x) = {V - AY{x) to all orders 
in the external momenta q flowing through the pseudoscalar current sources and 
Px = ^ + Q through the left current sources. This generalized four-point function 
will afterwards be reduced to a /S.S = 2 two-point function by integrating in 
the X-boson momentum with an explicit cut-off /x. We have to calculate for 
instance the generalized four-point function {Piis{0)A'^^{x)Pds{y)A.lf{z)) (also the 
one exchanging by V'^ sources). For definiteness, let us give one example 
for each of the two general type of contributions to this generalized four-point 
function. They correspond to the Figures in ^ and ^. In each of the four- 
fermion vertices any Dirac structure conserving the strong interaction symmetries 
is allowed. 

The three-point-functions like contributions (Fig. |^) to this generalized four 
point function {Pds{0)A'^^{x)Pds{y)Alf{z)) consists, then, of two full three-point 
functions (all orders in external momenta and quark masses) with one pseu- 
doscalar current source Pds and one axial-vector current source A'^^ each. They 
are obtained by gluing to the one-loop three-point functions full two-point func- 
tion legs (any permitted by the strong interaction symmetries) to obtain the full 
structure (see [0). Then the third leg of both full three-point functions is re- 
moved and the two three-point functions are pasted together with a propagator, 
i.e. a full two-point function with flavour either dd or ss plus a pointlike coupling. 
This propagator can have any Dirac structure compatible with the strong interac- 
tion symmetries. The flavours are also the ones corresponding to the generalized 
four-point obtained after conserving flavour in each four-fermion vertex. As an 
example, one contribution of this type to the generalized four-point function with 
flavour structure as indicated in Fig. R, is 



{PAPA) = 



g^p - gvU^pi-q - r 



ds ' 



X 



^PAv) {q + r,-r) Ig^c - gy'^-l^i-r] 



dd 



\ sdd r . 

^PAv) {-q-r,r) [g^s - gv'^^usiQ + r] 
+ ••• . 



ds 



(A.2) 
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The two- and three-point functions here are defined with the notation used in 
Ref. The non-barred n-point functions correspond to the full functions 



(all orders in external momenta and quark masses) and the barred ones to the 
one-loop expressions p4|, There are more than 320 contributions like this 
one. In these three-point-like function contributions one also can have products 
of two anomalous three-point functions, i.e. three-point functions which are pro- 
portional to a Levi-Civita symbol. To include them consistently we followed the 
prescription given in Ref. and in particular here we do not need to add any 
counterterm to the naive Feynman diagram calculation. 

The four-point like functions contribution to the generalized four-point func- 
tion consists, then, in full four-point functions with the same fiavour and Dirac 
structure as the generalized four-point function {Pds{0)L'^^{x)Pds{y)L^ff{z)). Each 
of these full-four functions is constructed by gluing to the one-loop four-point 
function two-pseudoscalar current sources Pds and two left current source L^f 
with the full two-point functions permitted by the symmetries of the strong in- 
teractions that gives the required structure. As an example, one full four-point 
function contribution to the generalized four-point function being considered and 
with the fiavour structure corresponding to Fig. ^ is 



{PAPA) = 



+ -[1+ gsllp{q)sd] 

, sdsd 



2 ' ^"'""^ Jo (27r)4 

X y^PAPA) {q + r,q,-q-r) \^g^s - gv^^siq + r 
+ ■■■ . 



ds ' 



ds 



(A.3) 



Here, the four-point functions notation follows up the notation of the two- and 
three-point function notation explained before and introduced in [p4| , [27[] in an 
obvious manner. There are 16 contributions of this kind. 

The one-loop n-point functions are regularized using a proper time cut-off 
that introduces the scale A^, see Refs. for details. Then, these two 



examples give idea of the kind of calculations one has to perform. We have 
performed several checks on our calculation. All of these checks of our ENJL 
model calculation, namely Ward identities, various comparisons with lowest order 
CHPT calculations,. . . , have been passed successfully, see Sects. ^, | and ^. 
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